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Faceted particles formed by the frustrated packing of anisotropic colloids on curved 

surfaces 
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We use computer simulations and simple theoretical models to analyze the morphologies that 
result when rod-like particles end-attach onto a curved surface, creating a finite-thickness monolayer 
aligned with the surface normal. This geometry leads to two forms of frustration, one associated 
with the incompatibility of hexagonal order on surfaces with Gaussian curvature, and the second 
reflecting the deformation of a layer with finite thickness on a surface with non-zero mean curvature. 

We show that the latter effect leads to a faceting mechanism. Above threshold values of the inter¬ 
particle attraction strength and surface mean curvature, the adsorbed layer undergoes a transition 
from orientational disorder to an ordered state that is demarcated by reproducible patterns of line 
defects. The number of facets is controlled by the competition between line defect energy and 
intra-facet strain. Tuning control parameters thus leads to a rich variety of morphologies, including 
icosahedral particles and irregular polyhedra. In addition to suggesting a new strategy for the 
synthesis of aspherical particles with tunable symmetries, our results may shed light on recent 
experiments in which rod-like HIV GAG proteins assemble around nanoscale particles. 


I. INTRODUCTION 

The formation of crystalline or liquid crystalline or¬ 
der on curved surfaces is frustrated, resulting in topo¬ 
logical defects. For example, while the densest packing 
of disks is a six-fold coordinated hexagonal lattice which 
perfectly tiles the plane, placing such a lattice on a spher¬ 
ical surface requires at minimum 12 five-fold coordinated 
disks [1, 2]. Similarly, placing a nematic phase of rodlike 
molecules on a spherical surface requires defects with net 
topological charge +2 [3, 4]. Despite intense theoretical 
and experimental research into ordering on curved sur¬ 
faces for more than 100 years {e.g. [1, 2, 4-28]), studies 
have focused on systems in which molecules are confined 
within the local tangent surface. 

Here, we study the coupling between particle shape 
and curvature that arises when anisotropic (rod-like) par¬ 
ticles tend to align along the surface normal. Our model 
is motivated by recent experiments on the assembly of 
rod-like HIV GAG proteins assembling around spher¬ 
ical nanoparticles [29], and the discovery that rodlike 
molecules with interparticle attractions can form stable 
monolayers [30-32]. We find that this geometry leads to 
two forms of frustration, one familiar from the packing 
of disks on a sphere, and the second reflecting the strain 
required to deform a layer of finite thickness on a curved 
surface. Frustrated sphere packing is an effect of intrin¬ 
sic geometry (metric in origin) and thus requires non¬ 
zero Gaussian curvature. In contrast, the strain associ¬ 
ated with layer thickness arises due to coupling between 
bend and splay and depends on extrinsic (mean) curva¬ 
ture. Thus, the latter effect occurs even in cylindrical 
geometries with zero Gaussian curvature. The interplay 
between frustrations of extrinsic and intrinsic origin leads 
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to a rich variety of faceted shapes, including those with 
icosahedral symmetry and other polyhedra with lower 
symmetry, whose features can be readily tuned by con¬ 
trolling template size and inter-particle interactions. In 
addition to elucidating how multiple forms of packing 
constraints can conspire to modulate order, our findings 
identify a new design strategy for synthesizing particles 
with aspherical shapes. Such particles are of great inter¬ 
est in the nanomaterials and colloidal communities, since 
they can assemble into an extraordinary range of struc¬ 
tures (e.g. [33-43]). Previous studies have shown that 
flexible ligands on isolated or clustered nanoparticles can 
order to form spherically asymmetric coatings [35, 44- 
75]. The current work demonstrates that rigid ligands 
can also order to form diverse morphologies that, despite 
their complexity, can be predicted from simple statistical 
mechanics arguments. 






Model I Model II 


FIG. 1. Schematic of the models. (A) Rods are modeled as 
rigid strings of beads. Examples are shown with nb = 6 beads. 
The bottom bead in each rod is attracted to the nanoparticle. 
The distance between the centers of the nanoparticle and the 
first bead of an adsorbed rod is denoted by Rnp. (B) Two 
inter-rod potentials are considered. In Model I, bead-bead 
attractions are distributed uniformly along the rods; in Model 
H, interactions occur only between the middle two beads, with 
the bead-bead interaction strength e increased to maintain the 
same net rod-rod interaction strength. 
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II. MODEL AND METHODS 


Spherical nanoparticle in 3D. Our system contains 
a rigid nanosphere with radius Rnp — crb/2 and A^rods 
rigid rods comprised of nb beads separated by a dis¬ 
tance equal to their diameter ab- Rods dynamically ad¬ 
sorb end-first onto the nanosphere, and are driven by 
lateral rod-rod interactions to form a layer aligned with 
the local nanosphere normal(Fig. 1). Beads from differ¬ 
ent rods experience excluded volume interactions repre¬ 
sented by a repulsive Lennard-Jones potential RLj(r; db), 
Eq. 1. In addition, beads with the same index (counting 
from n = 0 at the rod bottom) in nearby rods experience 
attractive interactions represented by a Morse potential 
Vmorse(?’;en,crb,a = 7.5/(Tb), Eq. 2, with e„ as the well- 
depth for the attraction between beads with index n: 
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where H{x) is the Heaviside function and Eghift ~ 
—0.001. We consider two distributions of attractions 
along the rods (Fig. 1): (Model I) uniform attractions 
with Cn = e Vn , and (Model II) attractions occur for 
the middle 1/3 of beads, with strength = 3e. All 
beads experience excluded volume interactions with the 
nanoparticle given by Vlj(^; Rnp) with R^p the center-to- 
center distance between the nanosphere and a bead on its 
surface. In addition, the first bead in each rod (n = 0) is 
attracted to the nanoparticle by Wiorse(^; Rnp, ck = 
7. 5 /( 75 ) with Cnp = 20 /cbT parameterizing the bead- 
nanoparticle interaction strength. 

Simulations. We performed Brownian dynamics sim¬ 
ulations with HOOMD [76-78]. The simulations used a 
set of non-dimensional units, described in Ref. [79]. In 
this article, we present all energies and lengths in units 
of the thermal energy /cbT and the bead diameter (Jb 
respectively. 

Simulation initial conditions were representative of a 
rapid quench from high temperature. Such a configura¬ 
tion was obtained in two steps. First, rod configurations 
and orientations were chosen randomly except that ex¬ 
cluded volume overlaps between pairs of rods or a rod and 
the nanoparticle were not allowed. Then, dynamics were 
integrated for ^ 5 x 10^ time steps with attractive inter¬ 
actions disabled (e = e^p = 0), but repulsive interactions 
operational. Once this initial condition was generated, 
attractive interactions were turned on, and the simula¬ 
tion was integrated for at least 2 x 10^ time steps, with 
time step size 10“^ dimensionless time units. 

Except where mentioned otherwise, all observables 
were measured after simulations had equilibrated. Equi¬ 
libration was assessed by monitoring the number of ad¬ 
sorbed rods and the number of clusters (defined below) 
on the nanoparticle surface (see Fig. 2 for an example). 



FIG. 2. Time evolution of the number of adsorbed rods (main) 
and the number of faces (inset) for a typical simulation of the 
spherical system, with nanoparticle radius Rnp = 7.0, rod 
length nb = 6, and bead-bead interaction strength e = 1. 
Throughout the article, lengths are presented in units of the 
bead diameter (Jb and energies in units of /cbT. 


We set the total number of rods Atot in a given sys¬ 
tem so that at most (^ 55%) of rods adsorbed onto the 
nanoparticle at equilibrium (ranging from Atot = 25 for 
Rnp = 1.0 to Atot = 1080 for Rnp = 7.0). This en¬ 
sured that, although the simulations were performed in 
the NVT ensemble, the driving force for adsorption did 
not diminish significantly over the course of a simula¬ 
tion. The results presented were obtained by averaging 
over several independent simulations for each parameter 
value (^ 10 for Rnp < 4.5 and ^ 5 for Rnp > 5.0). Fewer 
simulations were performed for larger nanoparticles be¬ 
cause computation times become prohibitively expensive; 
both the number of simulated rods and the equilibration 
timescale increase with nanoparticle size. 

Cylindrical model in 2D. To consider a case in 
which only mean curvature is present, we also analyze a 
simplified two-dimensional (2D) system, where rods are 
end-attached onto the outside of a circular boundary (ie. 
a 2D slice of a cylindrical nanoparticle). For simplic¬ 
ity, we do not simulate the dynamics of rod adsorption 
onto the cylinder. Instead, we consider a grafting den¬ 
sity equal to close packing on the surface of the cylinder, 
Wods = tt/ sin“^((7b/2Rnp), so the position of each rod 
end bead is fixed but rods can rotate. To investigate 
the large attraction strength (low temperature) limit for 
this system, we performed Metropolis Monte Carlo sim¬ 
ulations [80] in which randomly chosen rods were rotated 
by a small angular displacement (up to 0.01 radians), and 
the well-depth was gradually increased to e = 30 /cbT. In 
the limit of hard spheres and a short-range attractive in¬ 
teraction, 0 (^ 00 , the energy minimizing configurations 
of this system can be calculated exactly. This result and 
typical configurations from the Monte Carlo simulations 
are discussed in section III. 


III. RESULTS 

Except where mentioned otherwise, we will focus on 
the 3D system. All of these simulations resulted in 
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rods strongly end-adsorbed to the nanoparticle to yield 
a dense layer. In the absence of lateral attractions be¬ 
tween rods (e = 0) the bottom beads (which form the in¬ 
nermost layer) aquire local hexagonal order punctuated 
by 12 five-fold defects arranged with icosahedral symme¬ 
try, as observed for the packing of spheres on a spherical 
surface [5, 10]. 

A disordered-to-ordered transition occurs 
above a critical inter-rod interaction strength. 

For weak inter-rod interactions (small e), rod orientations 
are roughly random, leading to disordered arrangements 
of beads in the layers above the surface. Above a 
threshold value e*(nb,i?np), we observe a transition to 
an ordered layer in which rods are approximately parallel 
to their nearest neighbors. The value of e* is determined 
by the competition between the orientational entropy 
of the rods in the disordered state and the attractive 
interactions of parallel rods (Fig. 3). We approximately 
calculate how this competition depends on R^p and n\y 
at the end of this section. However, we first examine the 
unusual configurations that arise once the system enters 
the ordered state. 

Due to the nanoparticle curvature, rods cannot all 
maintain parallel orientations and fixed interaction dis¬ 
tances with their neighbors. Thus, the system under¬ 
goes a spontaneous symmetry breaking in which the sur¬ 
face is demarcated by line defects (or fissures) which 
form a polyhedron that inscribes the nanoparticle sur¬ 
face. Within each ‘face’ of the polyhedron, rods are 
roughly parallel and arranged with hexagonal close pack¬ 
ing. Interestingly, the defect lines do not necessarily form 
a regular polyhedron. 

To examine how curvature affects the surface layer ge¬ 
ometry, we measured the number of faces as a function 
of nanoparticle radius. As shown in Fig. 4, the number 
of faces increases slowly with i?np- Faces were identified 
by clustering neighboring rods with parallel orientations. 
To avoid calculating transiently disordered rods or small 
clusters as separate faces, we required a face to exceed 
a minimum cluster size nmin- The minimum cluster size 
must necessarily depend on nanoparticle radius, since the 
maximum possible size of a face decreases with i?np- For 
each nanoparticle radius, we identified a local minimum 
in the cluster size distribution, above which there were 
well-formed faces and below which there were individual 
rods or small clusters. Based on this, rimin varied from 2 
for Rjjp = 1.0 to 5 for i?np > 3. 

As shown in Fig. 5, the formation of line defects in 
the outer layer affects the packing within the innermost 
layer. Each of the 12 fivefold disclinations of the inner 
layer tends to be located at a vertex where several outer- 
layer faces meet. However, since there are typically more 
than 12 such vertices, not every vertex is associated with 
a disclination. Moreover, as the nanoparticle radius in¬ 
creases, we see that the disclinations extend into scars 
which extend along the polyhedron edges. The number 
and sizes of these scars are roughly consistent with the 
case of spherical particles adsorbing on the surface of a 



Interaction Energy, e/k^T 



FIG. 3. (A) Number of faces Waces formed as a function 

of the bead-bead interaction strength e for rod aspect ratio 
nb=6 and spherical nanoparticle radius i?np=4. Snapshots of 
final simulation structures are shown for indicated values of 
e. To aid visibility, the outermost beads are colored in red 
whereas others are colored blue. We determined e* as the 
value of the interaction strength for which diVfaces/de is maxi¬ 
mum, calculated by fitting the results for each set of (rib, Rnp) 
independently to a spline. The resulting transition interaction 
strengths are compared to results from the theories described 
in the text as a function of (B) rih for Rnp = 4 and (C) Rnp 
for nb = 6. 


nanoparticle (in which case the number of excess dislo¬ 
cations per scar would go as ^ 2R^p/a\) for 2R^p/a\) > 5 
[5, 10]). However, we observe that the inner-layer de¬ 
fects in the rod system are more disordered than those 
observed in the case of spherical particles. 

Despite the evident coupling between the locations of 
the inner-layer disclinations and the fissures in the outer 
layer, we find that the Gaussian curvature responsible for 
the former is not required for the formation of the lat¬ 
ter. In particular, simulations on the 2D cylinder system 
(which has only mean curvature) exhibited analogous fis¬ 
sures in the outer layer. Example low energy configura¬ 
tions for several nanoparticle radii are shown in Eig. 6. 

















4 




FIG. 4. (A) Snapshots of faceted structures in the spherical 

system shown for indicated values of the nanoparticle radius 
with Ub = 6 and e = 1. (B) The number of faces is plotted 
as a function of the nanoparticle radius i^np from the simula¬ 
tions (points) and as a function of ^ for the linear tilt theory 
(Eqs. (7), black line) and the theory accounting for bead in¬ 
tercalation (Eq. (4), red line). Error bars are estimated from 
5-10 independent trials at each parameter set. 


Factors that determine the number of faces. It 

is well known that incompatibility between a material’s 
preferred local packing geometry and long-range order 
can lead to defect-stabilized phases {e.g. Frank-Kasper 
phases, blue phases, and twist-grain-boundary phases), 
in which finite-sized unfrustrated domains are separated 
by regions of broken short-range order (defects) [2]. In 
our simulations, we observe an analog of this behavior on 
finite nanoparticle surfaces. The resulting defect geome¬ 
tries can be understood from the following simple models. 

Bead Theory: We consider a face (patch of nearly par¬ 
allel, hexagonally ordered rods) whose orientations are 
aligned with the sphere normal at the patch center, which 
we define as the z-axis (Fig. 7). Away from the center 
of the face, the nanoparticle curvature forces rods out of 
register. For the rod excluded-volume geometry that we 
consider (a chain of spheres), this displacement can be 
accommodated by intercalation of neighboring rods until 
the patch extends along the surface to an angle bounded 
by 

9a « 7r/6 + 0.5crb/-R„p, (3) 

where the second factor approximately accounts for rod 



FIG. 5. Simulation snapshots illustrating the relationship be¬ 
tween inner-layer disclinations and scars (left), and outer- 
layer fissures (right) for two different nanoparticle sizes 
{Rrvp — 4.0 (top) and i^np = 7.0 (bottom)). The images 
on the left show only the inner-layer beads, while the im¬ 
ages on the right show all beads. Each rod is colored accord¬ 
ing to the coordination number of its inner-layer bead, with 
colors(coordination) given by: red(5), blue(6), cyan(7) and 
white(8). Goordination numbers were determined by Delau¬ 
nay triangulation [81, 82]. 

granularity. Growth of the patch beyond 0^ would lead 
to a large strain energy which increases with patch size. 
This is avoided by terminating the patch and commenc¬ 
ing a new patch with a different rod orientation. The two 
patches are thus separated by a line defect. The number 
of faces Waces can then be estimated from the fraction 
of the nanoparticle surface occupied by a single patch, 
giving 

^faces =2/(1 — cos6>d) for Spherical nanoparticles (4) 

^faces =7r/^d for Cylindrical nanoparticles. (5) 

These relationships qualitatively agree with the numbers 
of faces observed in simulations (Figs. 4 and 6). 

For the 2D cylinder system, the bead theory exactly 
calculates the energy minimizing configurations in the 
limit of hard spheres and a short-range attractive inter¬ 
action, a ^ oo. The number of faces is given by 

= ceiling (iV,ods/floor(2 + V,ods/6)) (6) 

where the ‘floor’ and ‘ceiling’ operands account for the 
fact that the number of domains and the number of rods 
within a domain must be integers. Due to these con¬ 
straints switches between neighboring values as 

the nanoparticle radius varies (Fig. 6) and the resulting 
dependence of on nanoparticle radius is shown in 

Fig. 6. We confirmed the prediction of Eq. (6) by com¬ 
parison with Monte Carlo simulations (see section II). 
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FIG. 6. The number of faces in a circular geometry is plotted 
as a function of i^np for results from exact energy minimizing 
configurations (points) and the bead intercalation theory (red 
dotted line), and as a function of '0 for the theory (Eqs. (8), 
line). Snapshots from Monte Carlo simulations at several val¬ 
ues of Rnp are shown on top, with rods colored blue and the 
nanoparticle in green. 



FIG. 7. Schematic of a “face”, dehned as a cluster in which the 
constituent rods are nearly parallel. The base of the cluster in 
contact with the nanoparticle dehnes a patch that is bounded 
by an angle Od- 


Linear Theory: While the model just described de¬ 
pends on the local packing geometry of our rods, the 
effect is general. To show this, we consider deposition 
on a spherical or cylindrical substrate of a film with suf¬ 
ficient thickness 6 such that shear or tilt deformations 
correspond to lower energy than bending. For simplicity, 
we consider a linear tilt modulus G (although our simu¬ 
lated rod monolayer has a highly nonlinear stress-strain 
relationship). 

We then consider a circular domain of such a mate¬ 
rial as defined above. As the domain extends across the 
surface by an angle 0^, the material undergoes a strain 
along z given by ^{0) = tan^. The domain boundary 


gives rise to an interfacial energy characterized by a line 
tension cr, which we expect to be proportional to the film 
thickness 6. The total energy density 7/tot is given by 
integrating the strain energy over the domain area, and 
summing over the number of domains Waces- 

47ri7^ 6G ^-2- [(cos 6>d + sec (9d - 2) + V' "sm6>dj 

( 7 ) 


with ^ = R^^5G/(t the dimensionless ratio between shear 
modulus and line tension. The term in parentheses on the 
right-hand side of Eq.(7) gives the strain energy, the next 
term is the interfacial energy, and we neglect contribu¬ 
tions from the 12 disclinations required by the spherical 
geometry [83]. Because the tilt modulus and interfacial 
energy are both dictated by the strength of molecular in¬ 
teractions, we expect the line tension to be proportional 
to the film thickness thus, ^ oc Rnp/c^h with (Jb the 
molecular spacing. A similar analysis for a cylinder with 
radius R^^ and length L yields 


UcyM 

2'kR^^L5G 


1 


[(tan^d - ^d) + V’ 


(8) 


The size and the number of domains 0^ and Waces are 
obtained by minimizing Eq.(7) or (8). As shown in Eig.4, 
the linear theories qualitatively describe the simulation 
results, except that the predicted number of faces mono- 
tonically increases with R^^. In the simulated systems, 
the bead intercalation cuts off Waces at about 15 or 6 for 
spheres or cylinders respectively, and thus the bead the¬ 
ory is more quantitatively accurate for large i7np- It is 
difficult to see this cut off in the sphere geometry because 
simulations became computationally intractable for large 
i7np, but the effect is clear in the cylindrical geometry. 
In both geometries the number of domains approaches 
2 as -0 approaches zero (corresponding to small particle 
size or small tilt modulus). In this limit, the material in 
each hemisphere aligns with the z-axis with an equatorial 
defect between the two domains, resembling structures 
observed in simulation of alkanes on gold nanoparticles 
[69, 70]. The similarity between the results on spherical 
and cylindrical substrates emphasizes that the frustra¬ 
tion associated with depositing a thick film on a curved 
substrate is correlated to mean curvature, in contrast to 
the frustration of packing discs on a spherical surface as¬ 
sociated with Guassian curvature. 

An upper critical nanoparticle size for faceting. 
Eor sufficiently large nanoparticles the system will en¬ 
ter the thin film limit, with bending deformations lower 
in energy than tilt. (While T/gphere and Ucy\ monotoni- 
cally increase with i7np, the bending energy of an elastic 
layer is independent of R^p in a spherical geometry and 
decreases as R^^ for a cylinder.) We thus anticipate a 
critical nanoparticle size above which faceting does not 
occur, given by i7*p ^ with a~^ the width of the 

interparticle potential. 

Intra-bundle twist leads to another mechanism 
of frustration: While the interactions along the rods are 
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uniform in the simulations described to this point (Model 
I, Fig. 1), simulations with interactions confined to the 
middle two beads (Model II) resulted in an approximately 
regular dodecahedron. The 12 faced morphology arises 
because the rods form twisted bundles, which presumably 
are allowed due to the increased orientational freedom 
available to rod lateral interactions when the attractions 
are limited to the two middle beads. Inter-filament twist 
provides another form of frustration which can limit the 
lateral growth of bundles [84-87], which likely couples 
with the curvature induced frustration to decrease the 
critical face size. 

Given the magnitude of structural difference which 
arises between a seemingly small change in the interac¬ 
tion potential, we decided to further test whether these 
are equilibrium configurations. We performed an addi¬ 
tional set of simulations in which the system was first 
equilibrated under Model II interactions, then switched 
to Model I interactions and re-equilibrated, and then fi¬ 
nally switched back to Model II interactions. Note that 
this simulation protocol is not meant to model a partic¬ 
ular experiment (in which it might be difficult to switch 
interactions in situ), but rather to test whether the simu¬ 
lation outcomes depend on initial condition. As shown in 
Fig. 8, the 12-face structure obtained from Model II inter¬ 
actions reconfigures to a 9-face structure upon switching 
to Model I interactions, and then returns to the 12-face 
structure upon switching back to Model II. Since these re¬ 
configurations require extensive morphological changes, 
these reversible changes suggest that the structures in¬ 
deed correspond to equilibrium morphologies. 



-f 

— model II 

I _^_ 1 _^_L 

— model 1 

_^ ^ _ 


qL[-^^^-1-,-L_ 

^0 10 20 30 

Time Step /10^ 


FIG. 8. Reversible transitions between morphologies: the 
number of faces is shown as a function of simulation time 
step, for a simulation starting with the Model II potential, 
switching to the Model I potential, and returning to the Model 
II potential. The time series demonstrates that the simulation 
outcomes are robust to changes in initial condition. 

Calculation of disorder-to-order transition en¬ 
ergy. Finally, we consider how the order/disorder tran¬ 


sition interaction strength e* depends on the size of the 
nanoparticle and the rod aspect ratio. The formation 
of order is determined by the competition between the 
orientational entropy of rods in the disordered state and 
the attractive interactions of rods in the polyhedral state. 
The orientational free energy of a rod in the disordered 
state is given by /free/^eF = log 2, where we neglect ex¬ 
cluded volume with neighboring rods, assume that the 
nanoparticle surface restricts rod orientations to a solid 
angle of 27r, and that rods have free orientations in the 
reference state. To calculate the free energy of rods 
within a polyhedral face, we first estimate the rod-rod 
binding free energy from the partition function for the 
bound state of two neighboring rods: 

/3/b(e,nb) = - log[grod(e,^^b)/47r] 

27r 

grod(e,nb) = -—3j-^exp[(nb - l)e/kBT] 

^b^^morse Z-^n=l ^ 

( 9 ) 

with Qrod the partition function for the interaction be¬ 
tween two neighboring rods integrated over angular fluc¬ 
tuations, k/norse = 2e(a^ as the second derivative of the 
Morse potential evaluated at its minimum, and we have 
expanded the interaction potential to second-order. No¬ 
tice that the relevant rod-rod interaction strength pa¬ 
rameter is (nb — l)e, since the interaction between the 
innermost beads of neighboring rods is independent of 
their orientations. Each interior rod within a polyhe¬ 
dral face interacts with six neighbors, whereas those 
on the perimeter have on average 3.5 neighbors. We 
therefore approximate the binding free energy of in¬ 
terior rods as /b(3e), accounting for double counting, 
and perimeter rods as /b(3.5e/2). For a face with 
^rods rods, the fraction of rods on the perimeter can 
be approximated as /perim = 2(7r/A^rodsPHCp)^/^, with 
^rods = 47ri?npPHCp/c^b^faces(^np), PHCP = 2/\/3 and 
Afaces(^np) the number of faces as a function of nanopar¬ 
ticle size, given by Eqs. (3) and (4). The free energy per 
rod within a face is then given by 

A/(e, nb, Rnp)/kBT = (1 - /perim) /b(3e, nb) 

+ /perim/b(l-75e, nb) (10) 

Einally, the value of the transition e* is determined by 
A/(e*, nb, i^np) = /free- The calculated transition val¬ 
ues are compared to simulation results in Eig. 3 (B, G) 
as functions of nb and i?np- We see that the theoreti¬ 
cal prediction agrees closely with the simulation results 
except at small i?np, where the continuum limit breaks 
down. 


IV. CONCLUSIONS 

In conclusion, our simulations and simple theoretical 
models demonstrate that the geometric frustration in¬ 
trinsic to assembling anisotropic particles on curved sur¬ 
faces leads to diverse equilibrium morphologies. While 
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the specific number of faces obtained at a given param¬ 
eter set depends on the details of the interparticle inter¬ 
action potential, a simple model based on a generic lin¬ 
ear stress-strain relationship predicts qualitatively sim¬ 
ilar behaviors. Thus, it should be possible to generate 
polyhedral morphologies using a wide range of particle 
shapes and interactions. Extending our model to describe 
other interparticle interactions or substrate geometries 
would enable designing alternative assembly morpholo¬ 
gies. Controlling the strength of interparticle interactions 
(e.g. by temperature or depletant concentration) enables 
reversible switching between ordered, faceted morpholo¬ 
gies and disordered structures. Moreover, by changing 
the form of inter-particle potentials, it is possible to ob¬ 
tain highly symmetric structures {e.g. a dodecahedron) 
or asymmetric structures {e.g. the 9-face polyhedron). 

Our model was particularly motivated by two experi¬ 
mental systems. First, Dogic and coworkers [30] showed 
that a suspension of rodlike viruses and non-adsorbing 
polymer form 2D colloidal membranes, comprised of one 
rod-length thick monolayers of aligned rods. Deposition 
of such a monolayer onto a colloidal particle could be 
driven by depletion interactions [88] or complementary 
functionalization of rod ends and nanoparticle surfaces. 
It would be interesting to observe the morphologies of 
such assemblages as a function of nanoparticle size, deple¬ 
tant concentration, and depletant size. Second, cryoEM 
micrographs of HIV GAG proteins assembled around 


functionalized nanoparticles exhibit extensive scars that 
resemble features of the line defects and disorder observed 
in our simulations (Fig. 5). However, we note that the 
GAG protein assembles with a preferred curvature in the 
absence of a spherical template, in contrast to the rod¬ 
like particles considered here. The physics described here 
will be most relevant to cases in which the nanoparticle 
radius is small in comparison to the cone angle. For cone 
angles which are commensurate with the nanoparticle ra¬ 
dius, the frustration resulting from mean curvature would 
be avoided, and we would expect similar morphologies to 
those observed for packing of discs on a sphere [89]. 
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